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Abstract 



X 

Vh ' Electron gas in a wire connected to two terminals with potential drop is studied with the 



Schwinger-Keldysh formalism. Recent studies, where the current is enforced to flow with a 
Lagrange-multiplier term, demonstrated that the current enhances the one-particle-correlation 
function. We report that in our model, such enhancement is not guaranteed to occur, but 
conditional both on the potential drop and the positions where we observe the correlation 
function. That is, under a certain condition, spatially modulated pattern is formed in the wire 
owing to the nonequilibrium current. 

PACS codes(keywords): 05.70.Ln (Nonequilibrium thermodynamics, irreversible pro- 
cesses), 72. (Electronic transport in condensed matter), 85.30.Vw (Low-dimensional quantum 
devices (quantum dots, quantum wires, etc.)) 



Recently, Antal et al. have reported that nonequiUbrium current enhances the one-particle- 
correlation function contrary to our naive expectation that nonequilibriuni transport would 
destroy any orders ||l|, g, g]. They utilized a trick with which the nonequilibrium flow is in- 
troduced in a frame of equilibrium statistical mechanics. That is, they introduced a Lagrange 
multiplier A, and investigated the (equilibrium) ground state of the Hamiltonian given by 7i— A J, 
where TL denotes a model Hamiltonian of a wire and J denotes the space integral of the cur- 
rent operator. Namely, the current flow is regarded as an environment variable, which is to 
be given ad hoc. As for an explicit example of 7i, they took up the transverse Ising model 
and the XY model in one dimension, and found that the above-mentioned behavior, namely, 
the transport-induced correlation enhancement, occurs in these models. Cardy extended their 
work for a nonintegrable case (a lattice scaler field theory), and confirmed their observation [Q. 
This enhancement due to nonequilibrium transport had been studied extensively in the field of 
classical statistical mechanics [^ ^ [T], |^, |9|, |10|. It should be noted, however, that the quantum 
version has an advantage that the (nonequilibrium) dynamics is governed by the Hamiltonian 
itself, whereas the classical dynamics depends upon the kinetic rules which we implement. 

In this letter, we investigate the spinless electron gas in a wire, which is connected to two 
leads with different chemical potentials. Hence, in our case, the current is driven to flow by 
the potential drop as in actual experimental situations. We stress that we have employed the 
Schwinger-Keldysh formalism pll 03 , that allows us to treat the transport phenomena driven 
far out of equilibrium. Thereby, we report that the enhancement due to nonequilibrium flow is 
not guaranteed to occur, but depends signiflcantly both on the potential drop and the positions 
where we observe the correlation. As in the previous studies — note that the transverse Ising 
model and the XY model are equivalent to /ree-spinless-fermion models — we supposed that 
the wire electrons are free (quadratic). We believe that the many-body-correlation effect would 
not change our essential conclusions. 



Our model Haniiltonian is given by, 
^ = E (^fe - ^) ^i^fe+E (^'^ + '4-) RiRk+Y.t{clc^+i+4+,c,)+T'{cici+clcc.)+T'{cl,CL+c[cc.' 

k ^ ^ k ^ ^ i 

(1) 

The operator i[, (i?{,) creates a conduction electron with wave number k in the left (right) 
lead. Ck denotes the dispersion relation of these leads, and /i = eV/2 {—eV/2) give the chemical 
potential of the left (right) lead. Therefore, the parameter eV gives the chemical-potential drop 
between the leads. The scheme how we describe nonequilibrium flow is explicated afterwards. 
The operator c^ {ca') is the annihilation operator at the left (right) terminal; namely, c^ = 
(l/ViV) X]fc -^fc ^^^ *-«' ~ (Vv^)Sfc-^*: with the number of conduction electron levels N. 
The operator cj creates an electron at site i in the wire. We assume that the wire is sufficiently 
long, and the lead electrons are injected at the intermediate positions of i = 1 and L with the 
coupling (transfer) amplitude T'. We could implement the open-boundary condition for our 
wire, for instance, with use of the Fabrizio-Gogolin Tomonaga-Luttinger- liquid formalism under 
the open-boundary condition |Q; it would just cause inessential complications. 

Properties of quantum wire, especially, those on the electron conductance, have been studied 
extensively so far. It seems, however, that these studies have been rather inclined to the vicinity 
of equilibrium that can be managed within the Kubo formula, and little attention has been 
paid to the situation driven far out of equilibrium. The Schwinger-Keldysh formalism |0, n2J 
circumvents that restriction. Their formalism is demonstrated for nonequilibrium steady flow 
through a wire by Caroli et al. explicitly |l4| . We follow this description in order to investigate 
the nonequilibrium flow. In their description, the nonequilibrium steady flow is realized through 
infinite-time evolution from such initial state where both leads are disconnected from the wire 
(T' — 0), and are staying at each ground state with different chemical potentials. With respect 
to such time-evolved nonequilibrium state | ), the following Green functions are defined; 

G\^{i) = -i0(t)({c,(t),c]}) (2) 

G%{t) = ie{-t){{c,{t),c]}) (3) 

F^At) = -i([c.(i),4]). (4) 



They satisfy the foUowing respective Dyson equations, 

G^M - g^oj)+gUoj)^kiG\^{Lo) (5) 

Gf,M = 5|,M+5ffcMSMGr,H (6) 

F,j{io) = f,,{Lu)+gl,{cu)j:„iFi^{uj) + f,kico)^kiGt^{to), (7) 

The small-letter ones are those for the (unperturbed) initial state, and the self energy is given 
by Sfei — T'{Ska^ii + SkiSai + Ska'SLi + ^kLSa'i) [Q ■ It is notable that the Dyson equation of 
Fij is rather involved, and as a matter of fact, nonequilibrium characteristics are included in 

Quantity to be investigated below is the one-particle-correlation function, which is given by 
Cij = {cjcj +c^Ci). This correlation function is the central concern of the previous papers |l|, y|, 
where it is represented in terms of the spin language; (SfS^). In one dimension (Tomonaga- 
Luttinger liquid) , (long-range asymptotic forms of) all the correlation functions are interlocked 
each other. Hence, one particular correlation function, for instance, the one-particle correlation, 
is sufficient to gain the essential physics. (For example, when the one-particle correlation 
is suppressed, the density-density correlation develops instead.) The one-particle-correlation 
function is expressed in terms of the Green function Fij , 

a, = ^^{F^,{t^Q)+F,.,{t^Qj) (8) 

^L{F,,{u;)+F,,{coj). (9) 

Therefore, we need to know Fij. By means of successive use of the above Dyson equations, 
Fij can be expressed in terms of the unperturbed Green functions. Through a straightforward 
linear algebra, we arrive at the following expression; 

F^3 = h + 9aT'F^, + glLT'F^,, + hiT'C^^ + UlT'C^.^ , (10) 

with 

+9aaT'^9\L{9a'c.'T'ifLj + /liT'G^, + fLLG^,,,) + U'c'T'Glj)} (11) 



/ {(1 - 9ac.T"g{^)a - f^>^,r'glL) ~ 9'»c.T"9{L9l.'.'T'^9li} (12) 

Fa'] = {(1 ~ 9aaT' 9ll){9a' a'T' {f Lj + IliT'C^j + JllT'C^,^) + fa'a'T'Glj) 

+9l'o.'T'^9li{9lc.T'{h] + hiT'Gl^ + /ilG^,^) + ^oT'C?^-)} (13) 

/ {(1 - 9'c.o.T'^9\i){l - ffL'a'r'^ffLL) - 9laT'^9\L9l'a'T'^9li] (14) 

G'aj = 9aaT Glj (15) 

Gq'j = 9a'a''F G%j (16) 

, ^ g?,(i - giLr^^gg^aO + g?Lr^^ga^a^gi, , . 

[l- 9t{r'^9%o){^- 91lT'^91,o.')- 91lT'^91'c.'91iT'^9%o. ^ ' 

a ^ g£j-(i - gn^'^gao) + g^.^'^ggogii . . 

^^ (i-g?iT'2.9sJ(i-52Lr'2gS'a')-giLT''gS'o'g2iT'2gSo' ^ ^ 

Therefore, the remaining task is to express the unperturbed green functions exphcitly. As 
is seen above, informations of the terminals are involved within the local Green functions at 

a and a' . Here, we assume that the spectral property at the terminals is of the Lorentz type; 
namely, 

gUa'a'i^) = i/iu^ + iw) (19) 

goa,o'a'('^) = l/{uJ-iW) (20) 

F„„(w) = -2iWsgn{uj-eV/2)/{uj^ + W^) (21) 

i^a'a'(a>) = -2iWsgn{LU + eV/2)/{uj^ + W^), (22) 

with the band width of the conduction-electron spectrum W. 

As for the wire Green functions, we used the following expressions; 

G-Joj) ^ ^e"'^l^-^le-l"l/'^= (23) 
■' Vp 



27'"ip i:^|i-j| 



GlAuj) = ^e'^i'~^'e-l"l/"" (24) 

■^ Vp 

F,,{u) = -^cos(c.(*-j)/^'F)sgn(c.)e-l-l/-^ (25) 

Vp 

with the density of states p = l/(27r). Let us mention some notions about these expressions: We 
derive these expressions for a one-dimensional free-electron gas, assuming that the one-particle 
dispersion is linear, and the Fermi points are shifted towards the origin; uj — ±vpk. (We will 



assume that the units of time and length are identical (isotropic); v-p = 1.) These prescriptions 
are common to those of the bosonization scheme, and thus the validity is guaranteed especially 
as far as the low-energy physics is concerned. Moreover, it is notable that we readily treat a 
many-body-correlated wire, just substituting the above expressions with those of Luther and 
Peschel |lq] . The high-energy-cut-off factor ujc stands for the band width of the wire electron, 
and is related to the lattice constant a; ujc = upfcc = vpir/a. In the Tomonaga-Luttinger liquid 
theory, final expressions are obtained through setting k^ —> oo, because one is motivated to 
know the long-wave-length physics, where the lattice constant is renormalized to zero. Here, 
we retain the cut-off factor, because our wire length is supposed to be finite, and the length is 
to be measured by the lattice constant. We set the lattice constant as the unit of length; a = 1. 

The above formulae complete our scheme to calculate the one-particle-correlation function 
(g). In order to do the numerical integration of eq. (|9|), we adopted the Romberg algorithm 
|16| ; we proceeded the Romberg iteration until the output result converges up to the tenth 
digit. All the data presented below are calculated for the condition of the lead-electron-band 
width W — 2, the coupling amplitude T' — 0.5 and the wire length L = 20. 

Let us turn to present the numerical results. In Fig. m, we plotted the correlation Ci^s 
(H) against the voltage drop eV. According to the preceding reports, the correlation should 
be enhanced by the nonequilibrium potential drop eV 0, 0, 0. We found that in fact, the 
correlation Cy is influenced by eV . However, the dependence is not monotonic, but oscillates 
with respect to eV. In other words, our result indicates that it is not guaranteed whether the 
correlation is stabilized or not, but conditional on eV. 

In order to clarify the parameter range where this non-equilibrium effect becomes significant, 
we had swept various parameter ranges. Thereby, we observed that the oscillation period is not 
influenced very much by W and T' , but is proportional to 1/L. Hence, it is suggested that the 
oscillation period is governed solely by the energy scale ~ ujc/L. Note that the length of the 
wire does contribute to the energy scale, and thus the existence of the terminals is crucial in 
the present phenomenon. 



Secondly, we report that the enhancement of Cy depends on the distance i — j as well. 
In Figs. H and ||, We plotted the correlation C4,i over various distances i at eV — 0.2 and 
eV = 1.6, respectively. At eV = 0.2 (Fig. ||), we see that in almost whole distance range, 
the correlation function decays monotonically, and for long-range distance (i — j ^ L), to our 
surprise, the sign of the correlation alternates. That is, the correlation function suffers long- 
wave-length modulation. At eV — 1.6 (Fig. H), the alternation period becomes shortened. We 
observed that for larger eV, in general, the correlation alternates more rapidly. We stress that 
the modulation period is not related to the (equilibrium) Fermi wave length like the Friedel 
oscillation, because in our formalism the Fermi points are shifted towards k ~ 0. Hence, we 
conclude that the nonequilibrium current does shift the Fermi points out of fc = effectively. 
The Fermi-point shift is found to take place actually in the aforementioned models in refs. 
||l|, g|. In Figs. H and H, we notice that the correlation of eV = 1.6 dominates that of eV — 0.2. 
One might think that the correlation is enhanced by the nonequilibrium driving force eV. Yet, 
we found that this enhancement is rather accidental, and the correlation amplitude becomes 
suppressed gradually by eV; remember the behavior shown in Fig. nl 

Next, we show the short-distance correlation Ci.i+4 for various positions at eV = 1 in Fig. 
0. We found that the correlation is spatially anti-symmetric with respect to the center of the 
terminals. Therefore, the correlation-function data presented so far change the signs, if we re- 
define Cij so as to become closer to the opposite terminal instead; C'ij <-> CL-i+i,L-j+i- This 
anti-symmetry may be plausible, if we notice the following symmetry of our model; i ^^ L — i + 1 
and eV ^^ —eV. Namely, owing to the anti-symmetry, the net enhancement of the correlation, 
which is integrated over the wire, is always vanishing. This feature is precisely due to the 
existence of terminals, which are not taken into account explicitly in the previous studies. 

Finally, In Fig. g[ we plotted the electronic current / against eV. The current is calculated 
by means of the following formula derived by Caroli et al. [|^ ; 

I = -2^ / 2:^G'|iGii {{g^a - gaa)fa'a' + faaiOa'a' " 9a'a')} ■ (26) 

We found that the dependence of the current on eV is fairly monotonic; that does not exhibit 



any oscillatory behaviors. Hence, it is shown that our observations presented above are not 
directly related to the current flow strength. 

To summarize, we have calculated the one-particle-correlation function in a wire (||) which 
is subjected to two biased terminals. Being contrastive to the preceding report where the 
nonequilibrium flow is introduced ad hoc, our results show that the enhancement of the corre- 
lation function due to the nonequilibrium flow is conditional both on the bias and the positions 
where we observe it. In particular, at a certain condition, the nonequilibrium flow gives rise 
to spatially modulated structures. These features are related closely to the existence of termi- 
nals, which are implicit in the previous Lagrange- multiplier formalism. The situation that we 
had studied in this letter is quite realistic, and it may be realized in nano-technology devises. 
However, because we are concerned in the longstanding statistical-mechanical problem whether 
non-equilibrium transport enhances correlation or not, we concentrated deeply on the one- 
particle correlation function, that may not be measurable directly in experiments. It remains 
for future study to calculate quantities directly accessible in experiments. 
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Figure 1: The one-particle-correlation function C4^s (§|) for various chemical-potential drop 
eV. 
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Figure 2: The one-particle-correlation function C4^i (|8|) for various distances i over the wire 
for eV = 0.2. 
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Figure 3: The one-particle-correlation function C4^i (|8|) for various distances i over the wire 
for eV = 1.6. 
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Figure 4: The one-particle-correlation function Ci,i+4 (H) at various positions i in the wire for 
eV — 1. We see that the spatial distribution of the correlation enhancement is anti-symmetric. 
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Figure 5: Electric current against the potential drop eV. 



